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Abstract. The quantum control landscape is defined as the functional that maps 
the control variables to the expectation value of an observable over the ensemble 
of quantum systems. Analyzing the topology of such landscapes is important for 
understanding the origins of the increasing number of laboratory successes in the 
optimal control of quantum processes. This paper proposes a simple scheme to compute 
the characteristics of the critical topology of the quantum ensemble control landscapes, 
showing that the set of disjoint critical submanifolds one-to-one corresponds to a finite 
number of contingency tables that solely depend on the degeneracy structure of the 
eigenvalues of the initial system density matrix and the observable whose expectation 
value is to be maximized. The landscape characteristics can be calculated as functions 
of the table entries, including the dimensions and the numbers of positive and negative 
eigenvalues of the Hessian quadratic form of each of the connected components of the 
critical submanifolds. Typical examples are given to illustrate the effectiveness of this 
method. 
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1. Introduction 

The control of physical and chemical quantum mechanical processes has [1, 2] recently 
seen many laboratory successes, especially utilizing ultrafast shaped laser pulses as 
controls. In many experiments, adaptive learning algorithms are applied to seek optimal 
control fields [3, 4, 5]. These achievements collectively reveal a surprising ease in 
discovering effective control fields despite the presence of noise, imperfections and 
of severe constraints on the controls (i.e., limited bandwidth of ~20nm operating at 
~800nm central wavelength in virtually all of the experiments). This behavior suggests 
that some generic foundation must lie behind the ease of teaching a laser to guide 
quantum system dynamics to produce specified outcomes. Recently, the notion of 
quantum control landscapes was put forth to address this matter, where the landscape is 
the map from the space of control fields to some physical objective (e.g., quantum state 
transition probability [6, 7], the expectation value of a quantum observable [8, 9], or 
the fidelity of quantum gates [10]). The landscape critical topology (i.e., the topology 
of the set of critical points of the landscape) was analyzed to give insight into the 
effort to search for an optimal control. The initial exploration [6] of the landscape for 
the control of the statc-to-statc transition probability found that the landscape had 
only perfect extrema such that no false traps (i.e., local sub-optima in addition to the 
global optimum) exist for the control searches to be caught in, which provides a basic 
understanding of the ease of searching for optimal controls. 

The same feature exists for generahzed landscapes of the observable expectation 
value on quantum ensembles [8, 9] except for up to ~ A'"! saddle suboptima, where 
A'" is the dimension of the system. The details of these critical points depend on the 
observable and the system initial state. For real quantum systems, especially molecular 
systems whose electronic, vibrational and rotational states are involved, it is essential 
to investigate the affects of degeneracies upon the landscape topology to understand 
which classes of quantum control problems are expected to be successfully treated 
in the laboratory. Moreover, the landscape complexity needs to be considered when 
analyzing the scaling of control field search effort with N . This paper delves further into 
the control landscape topology for general finite-level quantum systems with arbitrary 
degeneracies in the statistical distribution of the initial state and the spectral structure 
of the observable. 

The paper is organized as follows. Section 2 presents a group theoretical analysis of 
the quantum ensemble landscape in terms of double cosets, as a generalization of torii 
obtained in [8]. Section 3 ascribes the determination of double cosets to the enumeration 
of contingency tables whose marginal constraints are the degeneracy indices of the 
observable and the system initial density matrix. Section 4 uses these contingency tables 
to calculate several important topological features of the landscape critical submanifolds 
and discusses the physical implications of the findings. Section 5 provides some simple 
illustrative examples, and conclusions are presented in Section 6. 
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2. The Quantum ensemble control landscape 

For a A^-level quantum system with non-dissipative dynamics, the evolution of its density 
matrix is described by pit) = U{t)pW (t), where p is the initial state and the unitary 
system propagator U{t) obeys the Schrodinger equation 

th^U{t) ^ H[e{t)]U{t), U{to)^I, (1) 

Here the evolution extends over the interval <t < T. The governing Hamiltonian H 
contains a time dependent control field e(t) that guides the evolution of the quantum 
system. The goal is to find control fields that maximize the expectation value of some 
desired quantum observable 9 of the system at the time T: 

J[e{-)]^tr{u{T)pU^{T)e}. (2) 

The optimization also can be defined for a non-Hermitian operator 9 [11], but here 
only the Hermitian cases will be studied. The quantum control landscape is defined as 
the mapping J: e{-) — > R from the set of all admissible controls into the real values 
of J. Accordingly, the optimal controls must be critical points of J in £(-)-space. A 
main objective of the landscape analysis is to identify the geometry of all possible critical 
points, as they specify the accessible stationary values of J upon variation of the control. 
A special issue of concern is the determination of whether false suboptimal critical point 
traps exist, as their presence could limit attempts to identify where to reach the absolute 
maximum of the landscape. 

As the landscape over the space of control fields is complex to analyze, a natural 
approach is to re-express the analysis on the unitary group: 

J{U) =tr{UpU^9), UeU{N), (3) 

whose critical topology is easier to determine. The relationship between (2) and (3) is 
manifested by the chain rule: 

5J^{VJ{U{T)),SU{T)), 

where VJ(f/(T)) is the gradient of J at U{T) in U{N), and 6U{T) is the variation of 
U{T) caused by a control variation Se{-) over [0,T]. The bilinear operation {A, B) = 
Retr(A^B) is a Riemannian metric defined on U{N). A control e{-) is called regular 
if any variation 5U (T) of the resulting U (T) is attainable by some admissible control 
variation, otherwise it is called singular [12]. For any regular control, the sufficient 
and necessary condition for 5J = upon any control variation is VJ(C/(T)) = 0, 
which implies that it is critical if and only if the corresponding U{T) is critical for (3). 
Moreover, the optimality status (i.e., as a minimum, maximum or a saddle point) of 
the regular critical points and their corresponding U{T) are identical (see proof in [13]). 
Singular controls can also become critical points of (2) when the corresponding gradient 
vector VJ{U{T)) is orthogonal to all attainable variations 6U{T), but in such cases 
VJ(t/(T)) is not necessarily vanishing. In this regard, the landscape critical points can 
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be classified as regular and singular critical points, respectively. For landscape studies, it 
is more important to assess the universal properties that are Hamiltonian-independent, 
which are essential for understanding the mounting successes in the control of various 
quantum systems. Hence we are mainly concerned with the set of regular controls, for 
which the reduction from £(•) to U (T) preserves the critical topology from the £(-)-space 
to U{N). The influence of singular controls will be contained in another work. 

The necessary and sufficient condition for [/ to be a critical point of (3) has been 
derived in [8, 14]. The basic concept is that, under the parametrization U{s, A) — e^^^U, 
where s e IR and = A, of a neighborhood of U in U{N), U is critical if and only if: 

^J[U{s,A)] =0, ^A^ = A, 
as s=o 

which gives 

iY{iA[e, UpU^]) = 0, V At = A, (4) 
and this is equivalent to 

[e, upu^ = 0, (5) 

which can be used to describe the critical topology [6, 8, 11]. For systems with non- 
degenerate p and 9. the critical submanifolds consist of A^! number of A^-torii embedded 
in the unitary group hl{N) [6]. The presence of degeneracies in p or ^ will merge these 
torii into a smaller number of larger critical submanifolds, and their explicit description 
will be resolved in the following sections. 



3. Chciracterization of the Critical Submanifolds 

Let R and S be the unitary transformations that diagonalize p and 9 respectively, i.e., 

p = R) pR = diag{\i, • • • , Ai; • • • ; A^, • • • , A^}, 

9 = S'^9S = diag{ei, • • • , ei; • • • ; e^, • • • , ^s}, 
where Ai > • • • > Ar are distinct eigenvalues of p with ni,---,nr multiplicities and 
ei > ■ ■ ■ > es are distinct eigenvalues of 9 with mi, • • • , multiplicities. The landscape 
functional can be written as 

J{U) = tT{URpR^U^S9S^) = tT[{S^UR)p{S^URy9] = tr{UpU^9) (6) 

where the automorphism U — S^UR also runs over U{N). Thus, without loss of 
generality, we can always assume that both p and 9 are diagonal. 

Theorem 1 Let U{n) be the product group U{ni) x • • • x lA{nr) where U{ni) is the rii- 
dimensional unitary group acting on the eigenspace of \, and lA{ra) — U{mi) x • • • x 
U{ms) is defined in the same manner relevant to 9. A unitary matrix U is a critical 
point of (6) if and only if it is in the double coset 

W(n)7rW(m) = {PttQ : P E U(m), Q e U(n)} 



of some permutation matrix tt. 
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Proof: Suppose that [/ is a critical point. According to (5), U must transform p 
into block-diagonal form, whose block lengths correspond to the degenerate subspace 
dimensions of 9. This block-diagonal matrix UpW can be subsequently diagonalized by 
a nil X ■ ■ ■ X nig block-diagonal unitary matrix P G W(m). Since unitary transformations 
will not alter the spectrum of p, the resulting diagonal matrix can be always written as 
TT^pn where tt is a permutation matrix that reorders the diagonal elements of p. These 
operations can be expressed as P^UpW P = irpn^ , which is equivalent to 

p = {7r^P^U)p{TT^P^U)l 

By definition, the matrix Q — tt^P^C/ must be an element of the stabihzer of p mU{N), 
which from group theory is identified with U{n). Thus, we obtain the decomposition 
U — PttQ. Conversely, one can easily verify that an arbitrary permutation matrix is 
critical, J{PnQ) — J(7r) and the critical condition (5) is satisfied for both tt and PttQ, 
implying that PttQ is also critical producing the same landscape value as tt. End of 
proof. 

Theorem 1 shows that the critical manifold is a union of double cosets: 

[j W(m)7rW(n), (7) 

tt£V{N) 

where V{N) denotes the permutation group over indices. The set Ai^^ = W(m)7rW(n) 
is the double coset of n in hl{N) with respect to tl{m) and lA{n). In the special case 
that both p and 9 are non-degenerate [8], the stabilizer of p and 9 are both products of 
N one-dimensional unitary groups, i.e., 

w(m) = w(n) = [war , 

the critical manifold as the double cosets can be verified to consist of A^! disjoint A- 
torii, each of which is labelled by a permutation matrix. Generally, the occurrence of 
degeneracies in p and 9 will "merge" two torii AItt and Ai^K' together if tt' happens to 
be in M.t^, thereby reducing the number of original disjoint critical submanifolds, but 
increasing their dimensions. These critical submanifolds can be taken as the generalized 
"torii". Let P be the set of all inequivalent partitions of V{N) with respect to the 
equivalence relation "~" on V{N). which is defined as vr' ~ tt if tt' G A4n- The 
characterization of the disjoint critical submanifolds can be identified as seeking all 
inequivalent double cosets of the permutation matrices in U{N) with respect to U{n) 
and U (m) 

M^\J U{m.)7cU{n). (8) 

TreP 

We now introduce the concept of a contingency table to simplify the abstract 
description of a critical submanifold expressed as the double coset of some permutation 
TT. Suppose that both p and 9 have their diagonal elements arranged in decreasing 
order. Let kij be the number of positions on the diagonal where the eigenvalues Aj and 
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ej appear simultaneously after imposing the permutation tt on 6. The contingency table 
consists of these nonnegative overlap numbers as arranged in Table 1, whose row sums 
are n^'s and column sums are m,'s. 



Table 1. Contingency table K. 





nil 


m2 


ms 


ni 


kii 


ki2 


kis 


"•2 


k2i 


k22 


kis 


Ur 


krl 


kr2 


krs 



For example, for a three-level quantum system where p = diag{OA, 0.3, 0.3} 
and 6 = diag{OA, OA, 0.2}, the corresponding marginal constraints are ni = 1 and 
n2 — 2 for row sums; mi = 2 and m2 — 1 for column sums. Suppose that a 
permutation transformation tt exchanges the second and third eigenvalues of 9, i.e., 
9' — Tr9n^ — diag{OA, 0.2, OA}. Then the overlap number of the first eigenvalue 0.4 of 
p with the first eigenvalue 0.4 of 9 is ku — 1, which appears at the first position on the 
diagonal of p and 9'; ^22 = 1 is the overlap number of the second eigenvalue 0.3 of p with 
the second eigenvalue 0.2 of 9, which can be seen from the second diagonal elements of 
p and 9'. Similarly, k2i = 1 and ki2 = 0. The resulting contingency table is shown in 
Table 2. 



Table 2. The contingency table for 
the illustrMti\ (' three-level system. 



1 1 

2 11 



Every permutation matrix leads to a contingency table, but the same contingency 
table may be produced from different permutation matrices which are mutually 
equivalent with respect to the relation "~" . Therefore, the contingency tables provide 
an equivalent description of the critical submanifolds: 

Theorem 2 Every critical submanifold of (3) is uniquely determined by a contingency 
table satisfying the above marginal conditions, and vice versa. The critical submanifold 
corresponding to a contingency table K can be expressed as the quotient set 

■^" = ""»fKr'' - where «(K) = n P) 

Proof: Define Ft^{P, Q) = PnQ, where (P, Q) G W(m) xW(n), as the bilateral operation 
on TT. The set of critical points that are equivalent with vr is characterized by the 
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stabilizer sta6(7r) of in W(m) xW(n), i.e., the set of matrix pairs {U,V) G W(m)xW(n) 
such that Ft^{U, V) = UnV = n. Hence, the critical submanifold A^^r can be identified 
as the quotient set of U{m) x U{n) divided by stablir), which is isomorphic to the set 
W(n) n TT'^U{ni)TT as follows 

stabin) = {(7^07^^ V) : V E U(n) n 7rtW(m)7r}. 

The intersection set W(n) fl 7r'''W(m)7r is a Lie subgroup of U{N) and can be 
decomposed into the product of smaller unitary groups as W(K) = U{kii) x ■ ■ ■ x U{kj.s), 
where the ki/s are the entries of the contingency table for tt. This leads to the expression 
(9). Therefore, every critical submanifold can be uniquely determined by a contingency 
table. On the other hand, given the contingency table shown in Table 1, one can shuffle 
the originally ordered diagonal elements of 9 such that kij of the rrij eigenvalues ej 
are moved to any of the Hj positions where the eigenvalues Aj of p are located. Any 
of such shuffles can be represented by a permutation operation corresponding to this 
table. Hence, any contingency table corresponds to at least one critical submanifold. 
In this manner, we prove the one-to-one correspondence between critical submanifolds 
and contingency tables. End of proof. 

The analysis above provides an easy means to analyze the critical topology 
of general quantum ensemble landscapes, by which the evaluation of the critical 
submanifolds on continuous Lie groups is reduced to a simple finite combinatorial 
problem solvable by a computer. Prom now on, we use the contingency tables to label 
the different branches of the critical submanifolds. The landscape value of the critical 
submanifold Mk corresponding to contingency table K is 

N r s 

1=1 1=1 j=l 

where tt is some permutation matrix whose contingency table is K. 
4. Chciracteristics of the critical submanifolds 

The contingency tables are not only convenient for labelling the critical submanifolds, 
but also powerful in identifying their intrinsic topological characteristics that are 
important for the performance of search algorithms seeking effective quantum optimal 
control field. These characteristics include (1) the dimension, Z)o(K), for a given 
critical submanifold Ma, which quahtatively reflects the size of the associated critical 
submanifold (especially the maximum of the landscape) that is crucial for assessing 
robustness to control field noise; (2) the numbers of the positive and negative "principal 
axis directions" near a critical submanifold, D_|_(K) and D_(K) [6, 8], which can affect 
the path taken by the search algorithm through the inffuence of the number of positive 
and negative directions in the vicinity of sub-optimal regions. 
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The dimension Dq{IC) in terms of the associated contingency table K can be 
calculated in an easier way by (9), with the fact that the dimension of U{t) is 

Do(K) = dim[W(n)] + dim[W(m)] - dim[W(K)] 

= + (10) 

1=1 j=l i=l j=l 

Computationally, the characteristics i^o(K) (resp., Z^+(K) and D_(K)) can be 
identified as the numbers of zero (resp., positive and negative) eigenvalues of the Hessian 
quadratic form (HQF) at any U e A^k, which determines the geometry in the vicinity 
of U. The HQF is defined as the second-order term of A in the Taylor expansion of 
J{e^^U) at [/, and it can be easily obtained as follows: 

n{A) = tr{AUpU^Ae - A^UpU^e). 

Let xp-y and yp-y be the real and imaginary parts of the /37-th matrix elements of A, 
which represent the coordinate variables in the tangent space of U. Then the Hessian 
form can be further expanded as the following sum [8]: 

= - E (^^-^7)(e^-e7)(47 + 2/|7)' (11) 

1</3<7<A^ 

where the real numbers Xg and {[i.^j = 1,---,N) are the rearranged diagonal 
elements of p and 9 after some permutation operation whose contingency table is K. 
Since there are always A*"^ terms in (11), the sum of the three indices should satisfy 

L>o(K) + D+(K) + L>_(K) = 

According to (11), the index is twice the number of (/3, 7) pairs for which 
{Xp — A^)(e/3 — e^) < 0, while D_ is twice the number of {(3,'~f) pairs for which 
{Xp — X^){efs — e-y) > 0. A critical submanifold is locally maximal (or minimal) if the 
corresponding index = (or D_ — 0), i.e., the Hessian eigenvalues are all negative 
(or positive). This happens only if the magnitudes of the A^'s and the e-y's are in the 
same (or opposite) order, which corresponds to a unique maximal (or minimal) critical 
submanifold. Therefore, there exist no local suboptima over the control landscape. 
Besides the absolute maximal and minimal critical submanifolds, all other critical points 
are saddles. In practice, the saddles will never form false traps for a search algorithm 
to approach global optimal controls, although the search may be slowed down when 
the algorithm runs in their neighborhoods. This behavior is a strong support for the 
observed relative ease of searching for optimal controls in quantum systems. 

Now suppose that Xp — Aj and A-y = Ap (1 < i,p < r); ep = ej and = 
(1 < J, ? < s)- By the nature of the contingency table, there are kijkpg possibilities for 
this to happen. Since Ai, • • • , A^ and ei, • • • , are both in decreasing order, we have the 
relationship: 

(A^ - X^){ep - e^) ^ ^ (i - p){j - g) ^ 0. 
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Hence, the D±{K) indices are equal to the following quadratic function of K: 



L'±(K) = 2 ^^^^^ 



(12) 



(i-p)0-9)^0 



which, by taking K as a r x s integer matrix, can be further written in a compact form: 



Here the entries of the index matrix = {cxij} (1 < i,j < t) are in upper triangular 
form 



We have not found an explicit formula for counting the number of critical 
submanifolds, which mainly counts the saddle critical submanifolds because there is 
always only one maximal and minimal critical submanifold. Generally, the breaking of 
the degeneracies in p and 9 (e.g., when p turns from a pure state to a mixed state) 
increases the number of saddle critical submanifolds, which may have adverse impacts 
on the search algorithms by creating a tortured path to approach the global optima. 
The examples in Section 5 analyze some special cases where the number of critical 
submanifolds can be explicitly calculated. For more complex situations, some good 
estimates have been found in [15, 16]. 

In real physical systems, there always exist disturbances or some other factors 
that may destroy the degenerate structures of the initial density matrix and observable 
operator. Since all the topological landscape properties are uniquely determined 
by their degeneracy degrees, even a very small perturbation may alter the critical 
topology. Heuristically, supposing the eigenvalues of p and 9 are well separated, the 
perturbation will wrinkle the existing "flat" critical submanifolds such that new smaller 
ones emerge, increasing the number of critical submanifolds up to A^!, while retaining 
only two absolute extrema. The submanifold dimensions and principal axis directions 
are consequently different in these cases. Nevertheless, provided that the perturbation is 
sufficiently small, e.g., compared with the step size taken upon numerical or experimental 
searching for a control, we expect that the structural changes of the critical topology will 
have little influence on the optimization algorithm performance, because the original 
critical regions will remain almost flat. To obtain a full view of the landscape, the 
geometrical curvature that reflects "steepness" (or "flatness" ) of the landscape near the 
critical submanifolds should also be considered. 

5. Examples 



L>+(K) = 2tr(J^KJ,Kt) 



(13) 

(14) 





In this section we will apply the results obtained above to several simple examples, to 
illustrate some basic features of the control landscape topology. 
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5.1. Pure state system 

Suppose that p is a pure state, then ni — 1 for Ai = 1 and n2 — N — 1 ior \2 — 0; the 
observable 9 has r distinct eigenvalues ei > • • • > with degeneracies being mi, ■ ■ • ,ms. 
The enumeration of the contingency tables is rather simple because the positions in 
only the first row can be filled with a single 1 and s — 1 zeros, which amounts to s 
different contingency tables as shown in Table 2. Therefore, there are s distinct critical 
submanifolds for such pure-state systems. Applying (10) and (12), the dimension and 
D±(Kj) indices of the j-th critical submanifold are 

DoiKj) = N{N -2) + 2mj, 

D+{Kj) =2(mi + --- + mj-_i), 

D^{Kj) = 2(mj+i + • • • + m,), 

for J = 1, • • • , s. One can specify that AIki is the maximum submanifold with 

Do(Ki) = N{N - 2) + 2mi, D-(Ki) = 2{N - mi); 

AIks is the minimum submanifold with 

Do(K,) = N{N - 2) + 2m„ D+(K,) = 2{N - m,). 

Obviously, high degeneracy in the largest eigenvalue of the observable may facilitate the 
optimal searches towards maximizing the cost functional and enhance the robustness of 
perfect control, while that of the smaller eigenvalues retards the control searching on 
level sets of lower landscape values. 



Table 3. Contingency 


tables K," 


for a pure 


state sys 


Umi. j = 1. ■ ■ ■ 








mi 




rus 


1 





1 





N-1 


mi 




nis 



5.2. Non- degenerate observable 

Suppose the spectrum of 9 is fully non-degenerate, i.e., nii = ■ ■ ■ — = 1; the density 
matrix p has r distinct eigenvalues Ai > • • • > with degeneracies being rii, - ■ ■ ,nr. 
Then, the entries in the contingency table can only be or 1, and there are 
choices for placing ni numbers 1 in the N positions in the first row. Having fixed 
the first row, there are — — rr choices to place Uo numbers 1 in the remaining 

' n2'-{N—m—n2y- f ^ o 

N — Ui positions in the second row, etc. Finally, the total number of distinct critical 
submanifolds is 

A^! (iV-m)! {N-m n^-2)! _ N\ 



ni\iN-ni)\ n2!(iV - m - 712)! 
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The corresponding dimension of each critical submanifold is 

r 
i=l 

Prom (9), these critical submanifolds are diffeomorphic to each other. This result 
generalizes the special case of both p and 6 being non-degenerate, whose critical 
submanifolds are all A^-torii. The D±{K) indices vary with the corresponding 
contingency tables, among which D_ = N"^ — at the maximum and = 

N"^ — at the minimum. 

5.3. A multi- degenerate system 

Consider an eight-level system with degeneracies n = (1,3,4) and m = (2,6). All of 
the possible contingency tables are listed as below: 








1 







1 




1 










3 


,K2 = 


1 


2 







3 




2 


2 




1 


3 




1 


3 




" 


1 " 




" 1 


" 








2 


1 


, K5 = 


1 


2 











4 







4 







The topological characteristics of the five critical submanifolds are summarized in Table 
4. The second critical submanifold is the largest with dimension 50. The minimum 
submanifold is the second largest. The maximum submanifold is the smallest, and from 
(9), it is diffeomorphic to the fourth saddle submanifold because they have the same 
group of nonzero entries in their contingency tables. The high degree of degeneracies 
yields the indicated high dimensions of the critical submanifolds. 

Table 4. Characteristics for the case of an eight-level system with degeneracies 



n= (1,3,4) and m= (2,6). 


No. 




1 


2 


3 


4 


5 


Manifold dimension 




48 


50 


46 


44 


44 


Positive axis direction 




16 


8 


6 


4 





Negative axis direction 







6 


12 


16 


20 


Type 




minimum 


saddle 


saddle 


saddle 


maximum 



5.4. Molecular systems 

A typical molecular system at room temperature and under strong field bound-state 
control usually possesses a very large number TV of accessible discrete levels. Suppose 
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the system has initially populated several non-degenerate vibrational levels l^;!), • • • , jv^), 
each of which has associated rotational states. Often the rotational states arc densely 
packed, and for illustration here, they will be taken as degenerate with rij {i = 1, • ■ ■ , r) 
associated with each vibrational state \vi). The population pi of \v,^ is then equally 
distributed over the rii rotational states. As the control field can import considerable 
energy into the molecule, N can be very large such that rij <C N . 

Consider the transition control from the initial state to a M-dimension subspace, 
where both M and N — M are larger than n^. The target observable is expressed as 
a projector 

1 ^ 
1=1 

on this M-dimensional subspace. The contingency tables corresponding to the critical 
submanifolds have to satisfy (r + 1) column-sum conditions rii, ■ ■ ■ ,nr, Nq, where 
Nq — N — Y7i=i ^wo row-sum conditions M and N — M, as shown in Table 

5. 



Table 5. Contingency tables for the 
molecular illustration. 





m 


rir 


No 


M 




kir 


kl^r+l 


N-M 


k2i 


k2r 





Since the contingency table is determined by the values of independent variables 
^iij • ■ ■ ) which can vary from to without violating the marginal condition, we 
can obtain the following number of critical submanifolds 

r 
i=l 

The global maximal submanifold corresponds to ku — rii, i — 1, ■ • • ,r, which from (10) 
has dimension 

Do^N^- 2(N - M)(ni + • • • + n^). 

Another special case is control to a specific final degenerate state. For simplicity, we 
assume that the degeneracies of the initially populated vibrational states are identical, 
i.e., rii — ■ • • — rir — n. The projector 6 then acts on a n dimensional subspace, 
i.e., M — n. In this case, the number of critical submanifolds equals the number of 
nonnegative unordered partitions of n in the first row of Table 4, which gives 

^ ^ (n + r)! 
n\r\ 

In this case, the dimension of the global maximal manifold is Dq = {N — -\- rt? . 
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The molecular illustration above has (i) a small number of initially populated states, 
(ii) a target population with a modest number of states and (iii) a very large number M 
of accessible states. These circumstances can arise quite commonly in molecular control, 
and they produce the outcome that the global maximum submanifold dimension Dq is 
very large and scales as ~ A^^. This behavior implies that finding a control on the 
maximum critical submanifold should be relatively easy due to the large size of the 
target submanifold and the lack of false traps. 

6. Discussion 

We have given a complete description of the critical topology of quantum ensemble 
landscapes for general finite-level systems over their propagator spaces. It is shown that 
the most important topological features can be calculated from a set of contingency 
tables that uniquely specify the critical submanifolds. This framework opens up the 
possibility to explore the landscape of more complex quantum systems. We have 
applied it to the control landscape of open quantum systems [13], and another potential 
application in the future is the analyses of control landscapes for an infinite dimensional 
quantum system as the limit of a series of finite dimensional systems, which can be dealt 
with using the techniques developed here. 

The landscape topology is based on a strong assumption that the system is fully 
controllable, which may be limited in practice. In particular, a broad class of partially 
controllable systems have their propagators constrained on a proper subgroup 5f of 
U{N), and the condition for some U e ^ to be critical is reduced to 

ii{iA[d, UpU^]) = 0, V A e g, (15) 

where g is the Lie algebra of $f . Obviously, every critical point of the landscape on 
IA{N) that belongs to still remains critical on ^ (there may exist additional critical 
points since the condition (15) is weaker than (4)). In certain cases, these critical points 
may degenerate from the original saddle points into local maximal (minimal) points. 
Hence, the results obtained here for fully controllable systems provide a basis for future 
landscapes studies of uncontrollable systems. 
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